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We study the collision of two slowly rotating, initially non boosted, black holes in the close 
limit. A "punctures" modification of the Bowen - York method is used to construct conformally 
flat initial data appropriate to the problem. We keep only the lowest nontrivial orders capable of 
giving rise to radiation of both gravitational energy and angular momentum. We show that even 
with these simplifications an extension to higher orders of the linear Regge-Wheeler-Zerilli black 
hole perturbation theory, is required to deal with the evolution equations of the leading contributing 
multipoles. This extension is derived, together with appropriate extensions of the Regge- Wheeler and 
Zerilli equations. The data is numerically evolved using these equations, to obtain the asymptotic 
gravitational wave forms and amplitudes. Expressions for the radiated gravitational energy and 
angular momentum are derived and used together with the results of the numerical evolution to 
q ■ provide quantitative expressions for the relative contribution of different terms, and their significance 

D ' is analyzed. 
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O ■ I. INTRODUCTION 

o 

' One of the most promising sources of gravitational radiation currently under investigation is the collision and 
coalescence of binary pairs of black holes of comparable mass [[p. The dynamics of such systems in the regime 
corresponding to large separation, where the interaction is weak, appears to be adequately described in terms of post- 

q-i Newtonian approximations At sufficiently close separation, where the interaction becomes strong, the system 
is expected to go through a "plunge" phase, eventually leading to the formation of a common horizon, after which 

(3jT)! the system behaves as a perturbed black hole that settles to a final stationary state. At this latest stage, "close 
approximation" perturbative method have proved capable of providing reliable information on waveforms and related 
. I quantities H. However, up to the present, no equally successful methods have been found that can handle the 
intermediate strong field regime, although some speculations have been presented, based mainly on extrapolations 
between the early and very late regimes. The fact that these methods may not reveal effects that may be unique to 
the strong field stage of binary interactions has been recently pointed out by Price and Whelan Q , in relation with 
the effect of tidal interactions in binary black hole inspiral, in the case where the black holes have moderately to large 
angular momentum. 

One of the particular close limit regimes considered in Q], that in which the black holes have equal, parallel spins, 
pointing in the same direction and perpendicular to the line joining the centers of the black holes, had not been 
analyzed in detail previously in the literature ||, and therefore only rough estimates are given there. In this paper we 
provide that analysis by considering the collision of two slowly rotating, initially non boosted, black holes in the close 
limit. The black holes have equal masses and equal parallel spins aligned perpendicular to the line joining initially 
the centers of the black holes. The expression "initially non boosted" refers to the manner in which the initial data 
is constructed. In our case we consider initial data obtained by the "puncture" || modification of the Bowen- York 
ansatz [Q. Mostly for simplicity, we include angular, but not linear, momentum in the conformal extrinsic curvature. 
In the close limit considered here, this represents a perturbed, non axially symmetric, state of a (slowly) rotating black 
hole, and therefore it is natural to apply the Regge-Wheeler-Zerilli ||, || black hole perturbation theory to analyze 
its evolution. However, as we show in this paper, even with the above simplifications the evolution equations for a 
part of the leading contributions are linear, but not homogenous, and require formulating the problem in second order 
perturbation theory JTo|, for a problem with more than one perturbation parameter. The corresponding evolution 
equations contain "source" terms, and their influence on the asymptotic behavior of the solutions, which contains 
the information on the gravitational waves, needs careful evaluation. Since these expressions had not been given in 
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the literature before, this paper contains the details of construction of the initial data and the evolution equations, 
together with the corresponding expressions for the gravitational wave amplitudes, and radiated energy, in terms of 
(gauge invariant) Regge- Wheeler- Zerilli functions. 

In the close limit the system analyzed here bears a marked resemblance with the "pseudoinspiral problem" considered 
in |Tl| , where the binary black hole system has orbital, but no spin angular momentum. This resemblance stems 
from the fact that in the close limit, once a single horizon has been formed, both systems behave as distorted black 
hole settling to a final Kerr form, with the distortion dominated by quadrupolar terms. Since the systems lack axial 
symmetry, it is expected in both cases that the gravitational waves will carry also radiation of angular momentum. 
This expectation was confirmed by the computations carried out in , [Q for the "pseudoinspiral case" . Since 
only linear perturbation theory is required in that case, the computations were done both in the Regge- Wheeler- 
Zerilli and the Teukolsky formalism |l3[ . However, it was rather surprisingly found that although the radiated energy 
computed in both formalism showed a reasonable agreement, this did not happen for the predicted radiation of angular 
momentum. In fact, except for very small values of the angular momentum, the Teukolsky formalism, as applied in 
, , leads to the rather disturbing result that the angular momentum increases as a result of the radiation (the 
angular momentum radiated has opposite sign to that of the system). This is not the case for the prediction given 
by the Regge- Wheeler formalism. Although the detailed reasons for this discrepancy are not clearly understood, it is 
the opinion of the present authors that it may be due mainly to the rather large departure of the initial data from 
that of a Kerr black hole. Since the system considered in this paper has initial data that is closer to that of a Kerr 
black hole, in the sense that the perturbations are of second order in the angular momentum parameter, rather than 
first order as in the case of |ll| , , one of the main purposes of our analysis was to compute the angular momentum 
radiated in the process, as the disturbed black hole settles to its final Kerr form, and its eventual comparison with the 
corresponding Teukolsky results. As is shown in this paper, within the Regge- Wheeler- Zerilli formalism, this, as well 
as the computation of the leading wave forms and radiated energy, requires an extension of the theory to second order, 
(in fact, to even higher order, if we derive the expression for the radiated angular momentum as indicated in Section 
VII), which is explicitly given here. It should be clear that at this point that the alternative treatment be based 
on the Teukolsky formalism, and, possibly, its extension to second order jl4|| , is a task requiring as much extension 
as the one reported here. For this reason, that development, and its comparison with the results obtained with the 
Regge- Wheeler formalism presented here, will be considered in a separate paper, currently under development. 

The paper is organized as follows. In the next Section we consider the construction of the initial data. This is 
naturally given in a "Misner" type gauge. However, as we indicate, this gauge is not appropriate for a perturbative 
computation, and we derive the gauge transformation to, and the form of the initial data in a "Boyer-Lindquist-Kerr" 
type of gauge for the £=1, odd parity perturbations, where these have a very simple coordinate dependence. We then 
consider the evolution of the second order perturbations in the Regge- Wheeler gauge, and obtain the corresponding 
Zerilli, and Regge- Wheeler functions and equations. We describe their asymptotic behavior, their relation to the 
corresponding gravitational wave amplitudes, and derive the expressions for the radiated power and energy. Finally, 
we obtain an expression for the radiated angular momentum directly using perturbation theory, in terms of the Zerilli 
functions for the perturbations. Comments and conclusions are given in the last Section. 



II. INITIAL DATA 



In applications of the Regge - Wheeler - Zerilli black hole perturbation theory one considers a family of metrics, 
depending on one or a few parameters, such that the Schwarzschild metric is recovered when the parameters are set 
equal to zero. The problem may be analyzed as an initial value problem, so that the specification of the members of 
the family amounts to identifying the corresponding initial data on a certain three dimensional hypersurface. Actually, 
since the computations are carried out in a coordinate patch that covers only the region outside the black hole horizon, 
one requires expansions of the initial that may be convergent only in that region. Examples are furnished by the 
families of initial data constructed via the "conformal approach" j7j] to the initial value problem in General Relativity. 
In this approach one assumes that the three-metric g a b is conformally flat, i.e. g a b — & 4 Vab, where $ is some non 
vanishing function, and r\ a b a flat three-metric, and defines a conformal extrinsic curvature K a b — § 2 K a i,. Then, 
assuming maximal slicing K ab g a b = 0, the initial value constraint (for a vacuum metric) are reduced to, 

V a K ab = (1) 
1 K ab K i. 

8 $7 w 
where all derivatives are with respect to the flat metric. In general, the conditions, 
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$ > , lim $ = 1. (3) 

p — >oo 

where p is a radial coordinate on the fiducial flat space where rj is defined, are imposed on <£>. Therefore, on the 
initial slice, asymptotically, for large p we have K ah — > iC' 1 '', and <? a (, — * rj a b- This fact is important because one can 
find exact solutions to the momentum constraint (Q), and, using (^), obtain a (partial) physical interpretation of the 
associated spacetime in terms of ADM observables, without solving (||) explicitly. For example, in M|, the solution 

~ 3 

Kab = — [tacdS c n d n b + n a e bc dS c n d ] , (4) 



where p = <J (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 , n a = x a / p, with x a cartesian coordinates on the flat space background, and S is 
a constant vector, is considered. Then, taking spherical polar coordinates for the flat background, so that 

r] ab dx a dx b = {dp) 2 + p 2 {def + p 2 sm(6) 2 (d(f>) 2 (5) 

one straightforwardly identifies S a with the angular momentum of the system, irrespective of the detailed knowledge 
of $. In the application given in | fL5f , the conditions (|^) are supplemented with appropriate boundary conditions for 
$, so that the resulting initial data is identified with that of a single spinning black hole (in a non stationary initial 
state), centered at the origin of the flat coordinate system. 

In this paper we consider the initial data that results from taking K ab as the sum of two terms of the form (Q), 
symmetrically centered around the coordinate origin, with the same S a . In the cartesian coordinate system associated 
with the flat background, we have, 

^ 3 

Kab = — [e a cdS c nfni b + n x a ^bcdS c n d ] 
Pi 

H — 3 [tacdS c n d n 2 b + n 2 a ^bcdS c n d ] (6) 
Pi 



where p x = ^(x 1 - L/2) 2 + (x 2 ) 2 + {x 3 ) 2 , p 2 = y/{x 1 +L/2) 2 + (x 2 ) 2 + (x 3 ) 2 , n\ = x a / Pl , and n a 2 = x a /p 2 . We 
also take S a = (0,0,5*). This may be thought of as the (conformal) extrinsic curvature corresponding to two black 
holes with equal angular momentum S, pointing along the z-axis, placed at the points (±L/2,0,0). Actually, this 
interpretation requires that <!> has a singularity structure in accordance with that of (|^). This may be achieved in 
different ways. Here we take the "punctures" ansatz Q, namely, we assume that $ is of the form, 

* = ®BL + $Reg (J) 

where, 

<S>BL = 1 + ^ + ^ 8 

2pi 2p 2 

that is, <&bl is taken as the Brill-Lindquist conformal factor ]lq ], and we impose that &R eg must be regular in the 
whole conformal plane, and vanish for large pi. In accordance with this prescription, <&R e g satisfies the equation, 

1 K ab K u 

with K ab given by (^|), which, for 5^0, may only be solved numerically. We notice, however, that for S = we 
should have $_R eg = 0, and, from the form ([)]), we expect that &R eg admits an expansion in powers of S 2 , near S = 0. 
Moreover, when L — > 0, &bl approaches the conformal factor for a Schwarzschild black hole. Therefore, the ansatz 
(0), together with (^), and (||), lead to a two-parameter family of initial data, with the Schwarzschild metric as the 
limit when the parameters vanish. In the slow, close approximation approach of this paper we consider the situation 
when the two holes are initially so close to each other that a single horizon is formed, and the system, from the point 
of view external to the horizon, may be considered as a perturbed single black hole. In this situation, the external 
field admits a multipolar expansion, and the multipolar order turns out to be correlated with different powers of the 
parameters, so that, by keeping the lowest orders in the parameters, we end up with a few multipolar terms, and this 
is expected to provide a good approximation to the exact initial data, and its evolution. 

There is always a certain degree of arbitrariness in the choice of perturbation parameters. Here we notice that, 
because of the conditions imposed on the initial data corresponds to a system with ADM angular momentum 
J = 2S. We shall therefore consider J, and L as expansion parameters. In the limit J = 0, L = 0, the conformal 
factor $ corresponds to that of a Schwarzschild black hole of mass M = 2Mq. In what follows we write all expressions 
in terms of M, rather than Mg. 
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A. Perturbative expansion of the initial data 



We recall that in the Regge- Wheeler black hole perturbation theory the analysis is restricted to the region outside 
the horizon. To construct our perturbative initial data we start by considering (|9|), but restricted to the region p > L. 
If we expand K a \, in powers of L/p, and keep only lowest orders, we find, 



K pp = -3 



JL 2 sin(20) sin 2 (0) 



p" 

9 JL 2 sin(2<^) sin(fl) cos(0) 
'8 7 



~ Jsin 2 (0) 3 JL 2 (7 + cos 2 (</>) (25 cos 2 (0) - 19)) sin 2 {0) 

p0 " 3 ^ 8 7 

K gg = ^JL 2 sin(2</>) cos 2 (0)/p 3 



Kb 



3 JL 2 sin(fl) cos(0) (1 + 3 cos 2 (0) - 5 cos 2 (0) cos 2 (0) 

4 ^ 

3 JL 2 sin(2<?!>) (4-9 cos 2 (0) + 5 cos 4 (0) 



= I — ( 10 ) 



This implies, keeping again only the lowest contributing order, 



K ab K ab = 18J 2 ^^ (11) 



Similarly, we expand &bl. 



M M(l~3cos 2 (0)sin 2 (0))£ 2 

$Si = 1+ 2^ ley (12) 



A simple power counting shows that the lowest order in $_R eff is J 2 , if we consider J 2 and L 2 of the same order. To 
this order (g) reduces to, 

The solution of this equation is, 

*fle S - $(0,0)(P)^>>^) + SfrOjO^M) (14) 

where Y"/™(0, (/>) is a standard spherical harmonic, and, 

, , 4V^J 2 (M 4 + 10pA/ 3 + 40p 2 A/ 2 +40 j o 3 A/ + 16p 4 ) _ C 2 
^ {P) 5M3(2p + M)5 + Cl + ~ 

, , \Z5^J 2 (M 4 + 10pA/ 3 + 40p 2 Af 2 +80p 3 M + 80p 4 ) „ 2 C 4 

(P) = 25p3 (2p + Af)5 + ^ + ? (15) 

We must set C\ = C3 = 0, to have the correct asymptotic behavior. The constants C2, and C4 are fixed so that 
^(o,o): an d ^(2,0) are regular for p = 0. This, for our perturbative expansion, is a consequence of the regularity 
conditions imposed on &R eg within the "punctures" ansatz || . See also the discussion given in |l7|] , for the analogous 
case of "boosted" black hole collisions. The final result is, 

40FJ 2 (Af 4 + lOpAf 3 + A0p 2 M 2 + A0p 3 M + 16p 4 ) 
™ (P) " 5A/ 3 (2p + A/)5 

_ . . 32^J 2 P 2 

^( 2 .o)(p)-- 25M(2p+M)5 (16) 
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This completes our perturbative computation of the conformal factor and extrinsic curvature. To construct the 
initial data for the evolution equations, we recall that if we fix the shift functions N a = 0, we have 

ds 2 = g ab dx a dx b - N 2 dt 2 
dg a b 



dt 



2NK ab (17) 



where for t = we identify g a b, and K a b with those obtained above. We next change from the conformal radial 
coordinate p to a Schwarzschild radial coordinate r, such that p = {\fr + \Jr — 2M) 2 /4, and choose N = yl — 2M/r, 
so that we recover the Schwarzschild metric, with mass M, for J = L = 0. The final results, in the standard Regge- 
Wheeler multipolar decomposition and notation M (see. e.g. JhJ, for more details) are listed below. Notice, however, 
that we indicate the odd parity perturbation components with k0i iTn , kli m and k2g iTn , instead of /i0^ m , hli t7n and 
/i2^ m , as in Q to avoid the explicit use of "odd" and "even" labels. 

B. Metric components 

The metric components obtained using the procedure just described may be written as follows, 

1. £ = 

4 y¥J 2 (rM + M 2 + 2r 2 ) 
H2 0i o = A , = - (18) 



2. 1 = 2 

H2 2 = K 20 = 4V ^ J2 16 (19) 
25 Mr 3 5 ^(^+Vr-2Af) 5 V ' 



™ _-^..- 8 ^ L2M ^ (20 ) 
K 2 K 22 8 ^ 2 fV30 (21) 



d t H2 2 _ 2 — 


—d t H2 2 2 - 


d t hl 2 _ 2 = 


8 iJ(8? 
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d t hl 22 = 


-d t hl 2 _ 2 


9tG 2 -2 = 


— 9tG 2: 2 = 


<9t-^2, -2 = 


-d t K 22 = 



64 iJy/^L 2 Vr- 2MV30 
~5~ r( 7 /2) (y7 + vV-2M)4 



r (5/2) (y9+ Vr-2M) 5 



16 iJ^/nL 2 Vr- 2MV30 
15 r (7/2) (y 7 ? + — 2 M) 4 " 

_16 iJy/^L 2 Vr-2MV30 
~ ~5~ r (7/2) ( y7 + vV - 2 M> 



(22) 

(23) 
(24) 
(25) 

(26) 



3. 1 = 1, odd 
dkl _ 4 V3J0r 16 V3J^L 2 M 

dtk l lt - 4 — — - y (Vr + Vr _ 2M)5r( 5/2) ^ 
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3, odd 



d t kl s _ 2 

9**13,2 

d t k2. 



4 J L 2 V^V210(25r + 25 VrVr -2M-6M) 
105 r (5/2) (^7? + Vr-2M) 5 

8 y^V7L 2 J(25r + 25 v / rv / r-2M-6M) 
"35 (V? + VF^TM) 5 r( 5 /2) 

4 JL 2 V^\/210(25r + 25 A A : \/7--2M-6M) 
105 



■ t fi,^ 3j _ 2 — ^t^3, 2 



r (5/2) + Vr- 2M) 5 
J^/nL 2 Vr- 2MV210 
21 r( 3 /2) Vr-2M) 4 

16 VWr -2ML 2 JV7 



dt.k2o , = — 



7 (y? + Vr-2M) 4 r( 3 /2) 



(28) 
(29) 
(30) 
(31) 
(32) 



C. Extrinsic curvature components 

The extrinsic curvature components may be written in the form, 



A-fcli.0 



^0 



J. £= l,odd 
16 

+ 



y/ZJypiiU-M 



r (3/2) v / r _2M 5 (y/F + - 2 Mf r 2 Vr-2M 



(33) 



iffl-22, -2 = 

Khl 2 ,- 2 = 
KG 2 . _ 2 = 
KK 2 . _ 2 = 



-KH2 2 , 2 = - 



2. 1 = 2, even 
32 iL 2 jV30y^ 



-Khl 2i2 
-KG %2 - 
-KK 2 . 2 - 



5 (Vr + Vr - 2M) 4 r 3 
4 i J (8 r + 8 Vr - 2 M - 3 M) L 2 V30 
15 r 2 Vr ~ 2M (Vr + ^JT~^2Mf 

8 iJ^L 2 Vm 



15 r 3 (Vr + Vr-2M) 4 
8 i J^L 2 V30 
5 r 3 (VF + V^^M) 4 



(34) 
(35) 
(36) 
(37) 



3. e = 3, odd 



Kkl 3 ,- 2 = Kkl 3t0 = - 



2 J (25 r + 25 Vr Vr - 2 M - 6 M) V^i 2 V210 



105 



Vr - 2 M (Vr + Vr - 2M) E 



__ 4 y/^V7L 2 J(25r + 25VrVr-2M-6M) 
3,0 _ 35 (VF + Vr - 2M) 5 r 2 V^IO? 



Xfc2 3i _ 2 = irt2 3;2 



4 JV^i 2 V210 
21 r (Vr + Vr — 2 M) A 



Kk2 3 = l ^ L2J ^ 



7 (Vr + Vr - 2 M) 4 r 



(38) 
(39) 
(40) 
(41) 



where we are using our modified Regge- Wheeler notation, and aJfat the beginning of each function indicates that 
we have a multipolar component of the extrinsic curvature, instead of one of the metric. 
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III. THE CHOICE OF GAUGE FOR THE 1 = 1, ODD PARITY, PERTURBATION 



If we consider the initial data for the I = 1, odd, perturbation, we find that if we set fcO^o 
setting the shifts to zero, the evolution equations have the solution, 



Mi o = 4 



0, consistent with 



(42) 



This in turn leads to solutions for the other amplitudes that grow in time. For our calculations it will be more 
convenient to choose a gauge where, 



VtFJ 
V3r 



(43) 
(44) 



The gauge transformation required to produce this change is first order in J, for the I = 1, odd parity terms, 
and introduces changes in the remaining terms of order J 2 , and JL 2 . A simple computation, using second order 
perturbation theory, as given e.g. in |l0| , shows that for the £ = 2, even parity, order J 2 terms, the initial data 
resulting from this change is the same as for the initial gauge except that now we have, 



HQ 2 ,o = - 



16V5vrJ 2 
15r 3 (r-2M) 



(45) 



Similarly, the transformation induces changes in the I — 2, even parity, order JL 2 terms. In this case, these 
correspond to first order perturbation, because the terms are linear in J. A straightforward computation shows that 
the expressions derived previously are replaced by, 



dtH2 2 2 
h0 2 ,2 

<9(/ll 2 ,2 



dtG 2 2 



-d t H2 2 _2 



16 i V^V30(V^ + 3\/^ 
~5 



2M) 



= -M) 2 



1G 



r (7/2) (^/f+ Vr-2M) 4 



15 r( 3 / 2 ) + Vr - 2 M) 5 



-d t hl 2 .-2 = 77" 



8 i^Vm 



-8 



2 M - 9 M) 



d t K 



t"2,2 



dtG 2 , _2 



-d t K 2 _ o = — 



15 yf + Vr-2Af) 5 r( 5 / 2 ) 

16 iVTrVMVr - 2M 
15 r( 7 /2) (V?+\/r-2Mj4 

16 i v 7 ^ n/SO (-4 M + V? - 2 M - 
5~ 



r (7/2) (^/f+ Vr-2Af) £ 



(46) 
(47) 
(48) 
(49) 
(50) 



where we now have a non vanishing shift term (h0 2j ± 2 ^ 0), and the data not explicitly indicated is zero. 

Finally, the gauge transformation equations for the t = 3, m = 0,±2, odd parity perturbations take the form, 



fcl 3 , 

fc0 3 ,±2 
&2 3 , ± 2 



/35r 
fcl 3 , - 

fc2 3:0 

2^7 r 

~21r 

kh,± 2 



K. 



2,0 



(it 



2G 



35r 2 



L 

K 2 ,± 2 ~ 
2t 



K. 



2.0 

L 2 
2,0 



L 

2^2,0 



2G L 



V7r 



2, ±2 



L 

^2,± 2 



2G L 



2, ±2 



= fc2 



3, ±2, 



(51) 



where K 2 ^ ±2 , G 2 ^ ±2 are the even parity, order L 2 , Regge - Wheeler functions in the Bowen - York gauge. In general, 
if necessary, in this paper we shall use a superindex L 2 , or L 2 J, to make explicit the order of the terms. The above 
equations may be used to give explicit expressions for the initial data in terms of r. We shall not display them here. 
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IV. THE EVOLUTION EQUATIONS FOR THE INITIAL DATA 



The time evolution of the initial data we have just constructed may be handled as follows. First we notice that 
the initial data depends on two parameters: J, and L. Actually, only even powers of L appear. This implies that, in 
a Regge- Wheeler type of expansion, terms of order J, or of order L 2 satisfy linear perturbation equation. We will, 
therefore, consider J, and L 2 as our perturbation parameters, and will restrict to the leading, non trivial, contributing 
orders and multipoles. In our case they may be classified as follows, a) Terms of order L 2 . These correspond to the 
I = 2, even parity, multipole, with m = 0,±2, and satisfy linear perturbation equations, b) Terms of order J that 
correspond to the £ — 1, odd, multipole with m — 0, and are related to the angular momentum. Their time evolution 
is known explicitly, c) Terms of order J 2 , with £ — 2, even parity, m — 0. These are second order in J, and obey linear 
equations with a "source" , and may be handled as in |l5| , but here we shall use second order perturbation theory to 
analyze their evolution, d) Terms of order JL 2 . These are actually higher order in perturbation theory. Since terms 
of order J have £ = 1, and are of odd parity, and terms of order L 2 have £ = 2, and are of even parity, in principle, 
they may contribute "sources" for terms of order JL 2 , and £ = 1,2,3, with both even and odd parity. However, an 
analysis of the form of the corresponding Einstein equations shows that the perturbations corresponding to I = 1,3, 
with even parity, and £ = 2 with odd parity, obey linear homogeneous equations. Since the initial data for these terms 
vanishes, they vanish for all times. Of the remaining terms, those with £ = 1 odd parity do not contribute to the 
radiation, and therefore, shall not be considered in what follows. In the next subsections we obtain and display the 
relevant equation for the evolution of the terms contributing to the radiation of angular momentum and energy, up to 
the orders considered. This implies, in some cases, that we shall have to extend previous results on the application of 
black hole perturbation theory in the Regge- Wheeler- Zerilli framework. These will be indicated where appropriate. 



A. Order L 2 



In this case we have to consider terms of even parity with I = 2, and m = 0, ±2. The reality of the metric implies 
that the amplitudes with m = —2 are the complex conjugate of the corresponding amplitude with m = +2, so we 
need to consider only one of these. Since all these perturbations obey first order equations, in the Regge- Wheeler 
gauge they may be uniquely written in terms of the Zerilli functions ip2,m(t, r), m — 0,±2, in the form fllOj, 

H0 2 , m (t, r)=H2 2 , m (t, r) 

d 2 ^ 2 ,m(t,r) , (2r 2 -2rM + 3M 2 )^ 2 , m (t, r) 



H2 %m {t, r) = (r-2M) 
-3 



dr 2 r(2r + 3M) 

(3 M 3 + 6 r M 2 + 4 r 2 M + 4 r 3 ) 



dr 



K 2 , m (t, r) 
Hl 2 , m (t, r) 



(2r + 3M) 2 r 2 
(6r 2 +6rM + 6M 2 ) 



r 2 (2r + 3M) 



i>2,m(t, r) + (1 



M, , d 



){g- r ^2, m {t, r)) 



(2r 2 -6rM -3M 2 ) 9^j, ro (t, r) .8 2 

+ r (^T^T V2,m{t, r)) 



(r - 2 M) (2 r + 3 M) dt v dt dr 

The Zerilli functions ip2,m(t, r) are given in a general gauge by the Moncrief fllq] expression, 



i>2,m(t, r) = 



r(r - 2M) 
3(2r + 3M) 

+ (2r + 3M) 



H2 2 , m -r- 



dK- 



2, m 



3M 



dr 



2M 



-K- 



2, m 



dG 



K 2 , m + (r- 2M) ( - ^hl 2 , m 

dr r z 



and obey the equation, 



d 2 i/j : 



where, 



'2, m _ d 2 ip 2 , m ij ~ 2M) (3M 3 + 6rM 2 + Ar 2 M + 4r 3 ) 
-g^— ~ dr *2 ~ 6 r 4 (2r + 3M) 2 ^ 2 ' m 

r* =r + 2Mln(r/(2M) - 1) 



(52) 



(53) 



(54) 



(55) 
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These functions are directly related to the gravitational wave amplitudes. The appropriate expressions are given 
below. 

The initial data are, 



i>2,o\t=o 



8V^r(7y/¥ + 5Vr - 2M)M 
l5(Vr + Vr - 2M) 5 (2r + 3M) 



dt 



= 



(56) 



and, 



^2,±2|t=0 



4V30^r(7Vr + Wr - 2M)M 
15(Vr + V7^2M) 5 (2r + 3M) 



2, ±2 



6>/ 







t=o 



and we notice that, since all the ip2,i(t,r) satisfy the same equation this implies 

[3 



ip2,±2(t,r) = 



4>2,o(t,r) 



(57) 



(58) 



To compute the wave forms, and radiated angular momentum and energy, we need to know the asymptotic behavior 
of the Regge- Wheeler gauge functions. This is immediately obtained from the knowledge of the asymptotic behavior 
of the corresponding Zerilli functions. Proceeding as indicated in |1Q), we obtain, 



^,m(*, r) = 4%(t r*) + 4%(t r*)/r + 41(t ~ r*)/r 2 + 4%(t - r*)/r 3 + G(l/r 4 ) 
where the functions ip-zlniz) ma Y be written uniquely in terms of a single function J-2, m in the form 



(59) 



2,m 



dx 2 



dx 



AW = 3^2,m - 3M 



dT- 



2,m 



-3MJT 2 , 



2\M 2 dF 2 ,r. 

4 da; 



(60) 



where m — 0, ±2, and the functions .7-2, m are determined by the initial data, through the evolution equations. 

The relation between the functions ip2,m and the gravitational wave amplitudes results from their gauge invariance. 
In an asymptotically flat gauge we find that, for large r, and t, we have 



i>2, m (t,r) ~ rG 2m (t,r) 



(61) 



where ~ implies equality up to terms of order r 1 . This result is used in Sections VI and VII, to compute the radiation 
of energy and angular momentum. 



B. Order J 2 



To order J 2 we only have perturbations for £ = 0, (that contribute only to the ADM mass), and with £ = 2, m = 0, 
of even parity, that play a role in the wave forms and radiated energy. The relevant part of the evolution of these terms 
may be obtained in different forms. One way is that considered in p^] . Here we describe an alternative treatment. 
As indicated previously, we fix the first order (order J) gauge, to the "Boyer - Lindquist - Kerr" form, klifl(t, r) = 0, 
kOifl = 4-\/ (n/3)J/r. This implies a second order (order J 2 ) gauge transformation on the original Bowen - York order 
J 2 data that, as indicated, introduces only a (time independent) change in the lapse function HO. 

Writing the second order Einstein equations it is easy to prove that the Zerilli - Moncrief function defined by, 



r{r - 2M) 

3(2r + 3M) 
^2 



H2o 



dK 2 



3M 



(2r + 3M) 



dr r~2M~ 
K 2 , + (r-2M) 



K. 



2. 



dG 2 



Or 



(62) 
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where H2, hi, K, and G are the corresponding Regge- Wheeler coefficients of order J 2 , obeys the equation, 



d 2 xlo _ d 2 xj _ (r - 2M)(3M 3 + 6rM 2 + 4r 2 M + 4r 3 ) j2 
dt 2 ~ dr* 2 r 4 (2r + 3M) 2 X2 '° 

16V5^(r - 2M)(9M 2 + l9rM + llr 2 ) 
5(2r + 3M) 2 r 6 



(63) 



where the second line contains the "source" terms, corresponding to the quadratic contributions of the order J 
perturbations. To establish the relation between this function and the gravitational wave amplitude we consider the 
asymptotic behavior of the perturbations for large t, and r. To this end we expand the perturbations in sums of terms 
of the form f n (t — r*)/r™, and obtain relations between the functions as in jl(|. The results, restricted to leading 
orders in 1/r, may be written in the form, 

J/2 i .„ = 12 G ° ( '~ r ' ) +0(l/r 4 ) 

M2i0 = _ 2 £VLzil +o(1/O 

G2 , = £V^ + 2 GW^) +o(1/r3) 

A - 2 , = 3 ^^ +6 evi_i! +o(1/rS) (64) 

where Gq(x) is a function of a single variable, determined by the initial data through the evolution equations, and a 
prime indicates a derivative. Replacing these expansions in (|62|), we find, 



xftfe r) = G" (t - r*) + ^ ^ + 0(l/r 2 ) 



(65) 



Therefore, asymptotically, to leading order in 1/r we have X2 o(^' r ) — r G'2,2^a(t,r), and xi ma Y be identified 
directly with the gravitational wave amplitude. The initial data for solving (^33|) is obtained from the results of Section 
II. After replacement and simplification, we find, 



j2 _ 4V57r(6r - 5M) 
X2,ol*=o — _ " 



75r 2 Af(2r + 3M) 

9x^0 



dt 

This is used in Section VI to obtain the indicated results 



= (66) 



C. Order J L 2 , 1 — 2, even parity, perturbations 

At order J L 2 we have I = 2, even, and £ — 3, odd parity, contributions. Formally, these are of second order, since 
we are considering both J, and L 2 as first order, and they satisfy linear inhomogeneous equations, with source terms 
bilinear in the order J, and order L 2 terms. Now, given any solution of these equations, together with those of order 
J, and L 2 , it is always possible to perform first a first order gauge transformation that takes the order J (I = 1, odd 
parity) terms to a standard form, that we have taken as the Boyer-Lindquist-Kerr form, given in Section III, and the 
order L 2 (£ = 2, even parity) terms, to the standard Regge- Wheeler gauge, where one may write the Regge- Wheeler 
gauge functions uniquely in terms of the Zerilli function, tp, (as in (|52|)), and the order L 2 Einstein equation are 
satisfied if ip satisfies the Zerilli equation. 

After the previous gauge transformations, a third gauge transformation may be used to put the £ = 2, even parity, 
order JL 2 terms, also in a standard Regge - Wheeler gauge form. In this gauge we have that one of the Einstein 
equations takes the form, 

8i dKkl 

m2 ^ m2 *-W^M)-ir w 
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where K is of order L 2 , and all other Regge- Wheeler coefficients are of order JL 2 . Taking this into account, it is 
straightforward to prove that, as in the cases considered in [flo|| , one can construct an infinite family of generalizations 
of the Zerilli function, that satisfy a Zerilli type inhomogeneous equation. From this family we have singled out the 
following function: 



r(r — 2M) r 2i 

X, 2 = 3^3^ [H2 2 , 2 - rd r K 2a] + -K 2>2 - ^ 2>2 (68) 



It satisfies the equation, 



d 2 X2 2 d 2 X 2 2 ( 2M\ (4r 3 + AMr 2 + 6rM 2 + 3M 3 ) 

6 1 - — Ia75Z I o»^2 ^2,2 



<9< 2 dr* 2 V r / r 3 (2r + 3M) 2 

8i(4r 3 + 56r 2 M + 36rM 2 + 15M 3 ) 9^2,2 
r 3 (2r + 3M) 3 dT~ 

The corresponding inversion formulas are, 



(69) 



(2r 2 - 2rM + Ml 2 ) 8x2,2 3(3M 3 + 6M 2 r + 4Mr 2 + 4r 3 ) 
r(2r + 3M) r 2 (2r + 3M) 2 



» ^ 9 2 y 2 2 (8r 4 - 60Mr 3 - 12M 2 r 2 + 15M 3 r + 27M 4 ) 9^2 2 

-l-O _ 9 Ml -i- 2?'- - — — 

1 ; 9r 2 r 2 (r-2M)(2r + 3M) 3 ) at 

„ (2r 2 + 16rM + 9M 2 ) d 2 ip 2 2 
+2i- 



3r(2r + 3M) 2 ) 9r<9t 

H§2,2 = H%2,2 — 



Si(r 2 + rM + M 2 ) dip 2 ,2 Ai d 2 ^ 2 



r 2 (r-2M)(2r + 3M) dt 3r <9r<9i 
= (2r 2 - 6rM - 3M 2 ) d X 2,2 d 2 X 2,2 
2,2 (r-2M)(2r + 3M) dt drdt 

2i(16r 5 - 192Mr 4 - 336A/ 2 r 3 - 564M 3 r 2 - 486Af 4 r - 135M 5 ; 



^2,2 



r 4 (2r + 3M) 4 
2i(16r 3 + 28Mr 2 + 78M 2 r + 45M 3 ) ^2,2 
+ 3r 3 (2r + 3M) 2 dV 

4i(r - 2M)(9M 2 + 18Mr + r 2 ) d 2 ip 2 , 2 
+ 3r 2 (2r + 3M) 2 dr 2 

_ 6(r 2 + rM + M 2 ) ( 2M\ d X 2,2 2i(2r 2 + WrM + 9M 2 ) dip 2 ,2 

2 ' 2 ~ r 2 {2r + 'iM) X2 ' 2 + ^ ~ r / + r 2 (2r + 3M) 2 dT ( ' 

Again, to compute the asymptotic behavior of the Regge- Wheeler gauge functions, we need to find the asymptotic 
behavior of the corresponding Zerilli functions. Proceeding again as in |10|, we obtain, 

X2,2(t, r) = X fUt - r*) + X ( 2 .kt - r*)/r + X { 2 l(t - r*)/r 2 + A - r*)/r 3 + G(l/r 4 ) (71) 



where now we have, 



A) d > ; 



dx 



2 



X2 ^ x >~ 6 dx 
(?\ , , dQn o d 2 T- 



2.2 



c?x da; 2 

\ IMC ■ 21M2 ^2,2 , .^2,2 , 23iM d 2 F 2<2 

and similar expressions for X2,— 2- The functions C/2,2 are uniquely determined by the initial data. This is obtained 
starting with the Bowen - York initial data, and first performing a gauge transformation of the order J, £ = 1, odd 
parity perturbations to a Boyer-Lindquist-Kerr gauge on the whole initial data, up to order JL 2 . This does not 
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modify the order L 2 terms but induces a (linear in J) change in the order JL 2 terms, and a second order gauge 
transformation on the order J 2 terms. On this modified data we perform a second gauge transformation, this time 
on the order L 2 , £ = 2, even parity perturbations to carry it to the Reggc - Wheeler gauge. This induces a further 
(linear) transformation on the order JL 2 perturbations. Actually, we do not need this last transformation because the 
homogeneous part of (|68|), (the part independent of ^2,2) is gauge invariant under gauge transformations of order JL 2 
of the £ — 2, even parity perturbations, and, therefore, we may obtain that part of the initial data without explicitly 
going to the Regge - Wheeler gauge. After some simplifications, the initial data for X2,2(i, r) may be written in form, 



X2,2|t=0 
9X2,2 







01 



t=o 



8iV30^(16r 2 + hlrM + 18M 2 - 15My/fy/r - 2M) 
15(Vr + Vr~2M) 5 (2r + 3M) 2 r 5 / 2 ' 



(73) 



We remark that the terms with I = 2, m — 0, even parity, of order J L 2 decouple from the order J perturbations, and 
satisfy homogeneous equations. Since the initial data for these terms vanish, they vanish and make no contributions 
to the evolution of the system. 

For the computation of the radiated energy and wave forms, we need the asymptotic behavior of the perturbations 
in an asymptotically flat gauge. It can be checked that in such a gauge, for large r, and t, we have 



X2,a(*,r) ~ rG 2 ,2(t,r) 



(74) 



where ~ implies equality up to terms of order r _1 , and therefore X2,2(i, r) is directly related to the gravitational wave 
amplitude. Similar results hold for X2,— a(*> >")• These are used in Sections XX and XX, to compute the radiation of 
energy and angular momentum. 

D. Order JL 2 , 1 = 3, m = ±2, odd parity, perturbations 

In this case we may consider again transformations that carry the metric to the Regge- Wheeler gauge form. The 
order JL 2 , £ = 3, odd parity, perturbations satisfy linear equations with a "source term" that is linear in the order 
L 2 , £ — 2, even parity, perturbations. In the Regge- Wheeler gauge these may be expressed entirely in terms of 
the corresponding Zcrilli function. One may in this case verify that if the Einstein equations are satisfied, then the 
Regge- Wheeler function defined by, 



Q 3 ,2(i,r) = -- 



2.. dk0 3 . 2 dkl 

-fc0 3i2 7T— "I 

r or at 



3.2 



105 



3^2 3 L , 

— h ~K 2 2 

or r 



(75) 



where K 22 is of order L 2 , and all other Regge- Wheeler coefficients are of order JL 2 , satisfies the equation, 



d 2 Q 3 ,2 d 2 Q 3 , 2 , 6(r-2M)(2r-M) 



dt 2 



dr* 



Q 



3,2 



where 



6(r 2 + rM + M 2 ) 
r 2 (2r + 3M) 



■lj'2 



4(r - 2M) 
V7r 4 

2M\ 



K L 



2.2 

r 1 dr 



(76) 



(77) 



Using these, and the Einstein equations, one finds the following "inversion" formulas, 



M) 3 , 2 = (r - 2M) 

r 2 



3Q 3 , 2 | (r-2M) 
dr r 

dQz.2 



Q3,2 + 



3\/7i 



■Ki 



(r - 2M) dt 

The asymptotic behavior of Q3,2 for large r, and t, is easily obtained from (|76|). It can written in the form, 



(78) 



Q 3 ,2 a Q {t - r*) + d(t - r*)^Q 2 (t - r*)l + 0(l/r 3 ) 



(79) 
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where 



Co(<) 



Qi(t) =6 



rf 2 Q(t) 
dt 2 
dQ(t) 



dt 



Q 2 (t) = 15Q(i) 



3MdQ(i) V7d 2 IC(t) 



dt 



21 a< 2 



(80) 



with Q a certain function that depends on the initial data, and JC related to V>2,2 by, d 2 K,{t — r*)/dt 2 ~ ip 2 ,2{t, r), for 
large r, and i. 

Using these results, and the gauge transformation that relates the Regge- Wheeler gauge to an asymptotically flat 
gauge we find that for large r, and t, we have, 



fc2 3 ,2(t, r) ~ -2rQ 3 , 2 (t, r) ~ -2rQ (i - r*) 



(81) 



The solution for m = —2 is the complex conjugate of that for m = 2. Therefore, solving (|76| we find the gravitational 
wave amplitude for this mode. 

The foregoing construction gives the complete solution of the perturbative Einstein equations for the mode consid- 
ered. For computational purposes, however, if one is interested only in finding the radiated energy and wave forms, it 
is simpler to notice that, after fixing the gauge for the order J ("Kerr" gauge), and for the order L 2 (Regge- Whheler 
gauge) perturbations, in a general gauge for the order JL 2 perturbations, the function, 



{3,2 



(r - 2M) 



fcl 



3,2 



r 2 d I k2. 



3.2 



(82) 



satisfies the inhomogeneous Regge- Wheeler equation, 

d 2 Q 3 ,2 d 2 Q 3 . 2 6(r - 2M)(2r - M) 



dt 2 



4(r - 2Mf 
V7r 5 



Q 3 ,2 



6(r 2 + rM + M 2 ) dip 2 ,2 



dtdr r(r-2M){2r + 3M) dt 



Asymptotically for large r, and t, Q 3 , 2 admits the expansion, 

Qs,a(*,r) = Qoit-n + Oir- 1 ) 



(83) 



(84) 



We remark that, after fixing the gauges for the lower order perturbations, Q 3 ,2 is gauge invariant under transfor- 
mations of order JL 2 . In particular, in an asymptotically flat gauge for these perturbations we have, 



and, therefore, 



fc2 3 , 2 (t, r), r ~ 2rQ 3 , 2 (i, r) ~ -2rQ (t - r*) 



k2 3 , 2 (t,r), t ~ -2rQ 3i2 (t,r) 



(85) 



(86) 



where ~ implies equality up to terms that decrease as 1/r for large r. This last equation establishes the relation 
between Q 3i 2 and the gravitational wave amplitude. 

To obtain the initial data for Q 3 ,2, we start with the Bowen-York initial data, and perform firt a first order 
transformation on the order J, i = 1, odd parity perturbations, that take these to the "Kerr" form, and then, on the 
result, we perform an order L 2 transformation on the i = 2, even parity, to carry these to the Regge- Wheeler gauge. 
Using the resulting form for the order JL 2 , I — 3, odd parity perturbations, we find that the initial data for Q 3 ,2 may 
be written in the form 



Q3,2(*,r)|t= =0 



dQ 



3.2 



dt 



4v / 2107i : M y/r-2M (5^? + lls/r - 2 M) 
105 (V? + \A- - 2 M) 5 rP 



(87) 



and we have the same initial data for m = — 2. 
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E. Order JL 2 , I = 3, m = 0, odd parity, perturbations 



In this case the Regge- Wheeler function is given by, 



i,o(*,i") = - 



10 



2 G>fc0 3 ,o , dfcl 3 ,o 

-fc0 3i0 7T— H ^r— 

r or at 



175 



dr 



-K. 



2.0 



(88) 



where again is of order L 2 , and all other Regge- Wheeler coefficients are of order JL 2 , and the Regge- Wheeler 
equation takes the form, 



d 2 Q 3 ,o _ d 2 Q 3 ,o 6(r-2M)(2r-M) 
<9t 2 dr* 2 r 4 



?3,0 + 



12(r-2M) A . £ 



35r 4 







where 



T ^ 



6(r + rM + M 2 ) 
r 2 (2r + 3M) 



2,0 



1 _ 2M\ chhfi 
r J dr 



and the "inversion" formulas are, 



WM = (r - 2 M)22M + <i^M) Oj , 



kh,o = 



dr 
dQ 3 ,o 



35r 



-Kk 



(r - 2M) dt 



As in the case m = ±2, the asymptotic behavior of Q 3 ,o for large r, and t, is easily obtained from 
written in the form, 

Q 3 ,o Co(< - r*) + Qi(t - r*)-Ga(* - + ^(l/r 3 ) 



(89) 



(90) 



(91) 
!)). It can 

(92) 



where 



Qi(t) = 6 



rf 2 Qft) 



dt 



Q 2 (t) = 15Q(i) 



3MdQ(t) 1 <9 2 /C(i) 



35 di 2 



(93) 

^2,0 



with Q again a certain function that depends on the initial data, and JC related to -02,0 by, d 2 K,(t — r*)/dt 2 
for large r, and t. 

Just as in the case m — ±2, we find that for large r, and t, the solution of (p9| ) is related to the amplitude in an 
asymptotically flat gauge by, 



fc2 3) o(*, r) =s -2rQ 3) o(*, r) =s -2rQ (i - r*) 



(94) 



Thus, again the solution of (|39|) provides the gravitational wave amplitude for this mode. 

As in the previous subsection, if one is interested only in finding the radiated energy and wave forms, we notice 
that, after fixing the gauge for the order J ("Kerr" gauge), and for the order L 2 (Regge- Wheeler gauge) perturbations, 
in a general gauge for the order JL 2 perturbations, the function, 



i,o(*,r) 



(r - 2M) 



kl 



3,0 



f fc2 3 , 
2 dr [ r 2 



(95) 



satisfies the equation, 



d 2 Q 3 ,o 5 2 Q 3 ,o 6(r-2M)(2r-A0 



dt 2 



dr* 2 

12 (r - 2M) 2 
35r 5 



Q 3 ,o 



9^2,0 , 6(r 2 + rM + M 2 ) ^2,0 



dtdr 



r(r - 2M)(2r + 3M) dt 



(96) 
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Asymptotically for large r, and t, Q3,o admits the expansion, 

Q 3l o(t,r) = Q {t-r*) + O{r- x ) 



(97) 



We also have in this case that, after fixing the gauges for the lower order perturbations, Q3,o is gauge invariant 
under transformations of order JL 2 . In an asymptotically flat gauge for these perturbations we have, 



and, therefore, 



fc2 8 ,o(*, r), r =s 2rQ 3 ,o(i, r) =s -2rQ (t - r*) 



fe2 3 , (t,r), t ~ -2rQ 3 ,o(*. r) 



(99) 



which provides the relation between Q3 ; o and the gravitational wave amplitude. 

Proceeding as in the case m = ±2, the initial data for Qs^ may be written in the form 



?3,o(*,r)|t=o = 



dQz,o 



Of 



t=o 



35 (Vr + Vr - 2 JV/) 5 



(100) 



V. EVOLUTION AND WAVEFORMS 



All the relevant Regge- Wheeler and Zerilli equations were evolved using straight forward modifications of the rout ines 
used in p9[ . As indicated there, the evolution is carried out in the Cauchy domain of some appropriately chosen interval 
of the r*-axis, and therefore, we do not need to impose any boundary conditions. To obtain a definite wave form we 
consider an extraction point at r of the order of 100M, (the exact point is actually irrelevant, since the amplitudes 
are independent of r, for large r), and register the different amplitudes as functions of integration time. For the £ — 2, 
orders J 2 and L 2 , even parity terms we obtain the familiar "quasinormal ringing" wave forms. This is also true for 
the order JL 2 , odd parity, £ = 3, amplitudes, although the ringdown frequencies are different from those for £ = 2. 
For the £ = 2, even parity perturbations of order JL 2 , we obtain a wave form of the same frequency as that for £ = 2 
of order L 2 , but with a different initial shape, and a noticeable shift in phase. This last feature is essential for the 
radiation of angular momentum. As illustrations of these results we include in Figure 1 a comparison of the wave 
forms discussed above. Notice that for \ we have indicated the imaginary part (the real part vanishes in this case) 




40 60 80 100 40 60 80 100 

t(M) t(M) 

FIG. 1. In the figure on the left the solid curve represents (in arbitrary units) the wave form corresponding to ^2,2, and 
the dashed curve that of Qs,2, as functions of time, for fixed r. Notice the difference in ring down frequency, for the different 
values of I. The figure on the right compares, again in arbitrary units, the wave forms corresponding to ip2,2 (solid curve), and 
to X2,2 (dashed curve). Here we notice the large difference in phase, giving rise to radiation of angular momentum. 
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VI. RADIATED ENERGY 



To compute the power and radiated energy we consider the metric written in a transverse traceless, asymptotically 
flat gauge. These conditions are satisfied asymptotically in what we have called a "Misner" gauge in this paper. We 
may therefore apply the formula |2C|] 



dPower . 1 
= hm 



16irr 2 



1 dh eq 
sin 9 dt 



1 fdh e 



4 a 



sin' 



9 dt 



Using the general form of the Regge - Wheeler expansion given in 0], for large r, we have, 



(101) 



188 



sin 2 i 



sm( 



r 2 E E 

1=1 m=-t 
oo m—i 

+ 2 E E k2 < 

1=2 m=-l 

oE E k2 ^ 



d 2 Yi 



1 d 2 Y, m cos 6> ay/ 



ae> 2 



sin 2 (9 



sin (9 36» 



1 <9 2 y m cos((9) 9F/ 



=2 m=-£ 
oc m—£ 

E E G ^ 

?=2 m=-i 



wn.(9)d6 6<t> sin 2 (9) 

a 2 Y f m cos 6» a y m a 2 y m 



sin 2 6»a0 2 sin 6» 36* 89 2 

d 2 Y f m cos 9 ay/ 



sin 000 , 



sin 



(102) 



where the sums over £, and m are restricted to the those of the perturbation considered. After expansion of the 
coefficients in the form, 



G t , m (t,r) = J 2 G^(t,r)+L 2 G^\t,r) + JL 2 G^\t,r) 
k2 t , m {t,r) = JL 2 k2 { l^\t,r) 
where we have explicitly indicated the orders of the terms, an integration over the angles gives, 



(103) 



3r 2 

Power = lim 

r-+oo Sn 



J 4 (gW 



2J 2 L 2 G ( 2 J ) G ( 2 L ) +LU [G 



2,0 



ZLt 2,2 ^2,-2 



5 J 2 L 4 / / . 



9 r2r4 ( A(Ji 2 ) ( A(Ji 2 ) 



(104) 



where an overdot indicates a/<9£. Therefore, in terms of the corresponding Regge - Wheeler and Zerilli functions, the 
total radiated energy is given by, 



Energy = lim — 
r— >oo Sir 



J 



(*g ) ) a+2jax9 ^a ,) +^ 4 («" ) ) s 



+ 20J 2 L 4 ( (Q^P) 2 + 2G,r ) + 2J 2 L^f ) X < 2 % ) 



dt 



(105) 



After numerical integration of the evolution equations, with the initial data derived previously, we obtain 

Energy/AT = 7.8 x 10~ 4 (J/Af 2 ) 4 - 2.8 x 10~ 5 ( J/M 2 f{L/M) 2 + 9.8 x 10~ 5 (L/M) 4 

+ 1.3 x 10~ 3 (J/ M 2 ) 2 (L/M) 4 (I = 2) + 2.1 x 1Q- 5 {J/M 2 ) 2 {L/M) 4 {£ = 3) 



(106) 



where in the second line we have written separately the I = 2, even parity, and i = 3, odd parity contributions of 
order JL 2 . 
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VII. RADIATION OF ANGULAR MOMENTUM 



The lowest order in the radiation of angular momentum comes from the interference of the £ = 2, m — ±2, even 
parity, outgoing waves of order L 2 , with those of order JL 2 , giving a contribution of order JL 4 . We may obtain the 
expression for the radiated angular momentum using the formulas in pjj . This implies, however, that we need to 
obtain the asymptotic transverse traceless gauge amplitudes, which implies a rather lengthy computation. We may 
instead use the fact that the following expression, 



J 



lim 

r — >oo 



1 



4V37T 



(r 2 fcOi, ,r (r, t) - 2rfc0i, (r, t)r 2 - kU,o,t (r, *)) 



(107) 



where fe0io(r, t), and fcli ) o(r, t) represent the corresponding Regge- Wheeler coefficients in an asymptotic expansion 
of the metric for large r, corresponds, for t = 0, and to leading order, to the (ADM) angular momentum on the initial 
hypersurface. It is clear that J, defined by ( |107| ) is time dependent. Since J{t = 0) represents the initial value of 
the angular momentum of the system, we should consider limj—joo Sf{t) 1 as the final value of the angular momentum, 
after the system has settled to its final Kerr black hole. The radiated angular momentum is therefore, 



AJ = J(0) - lim J(t) 



(108) 



To compute AJ perturbatively, we expand k0i^(r, t), kli^ir^t), in powers of the perturbation parameters J, and 



L. 



fc0i l0 (r, t) = J fc0$(r, t) + JL 2 fc0$(r, t) + JL 4 fc0$(r, t) + J 2 L 4 fc0^(r, t) + 
feli,o(r, t) = J kl{ 1] Q (r, t) + JL 2 kl[ 2) {r, t) + JL 4 kl[ 3) (r, t) + J 2 L 4 fcl^(r, *) + 



(109) 



Notice that there are no contributions of order J 2 , L 2 , and L 4 . If we consider now the Einstein equations for the 
different orders, we have that, to order J, 



J 



(i) 



■= (r 2 k0 {1 \ r (r,t) - 2rfc0 (1) (r,i) - r 2 kl {1 \ t {r,t}) 



4V3tt 



(110) 



is a gauge invariant constant, that, on account of the form of the initial data, is equal to J on the initial hypersurface. 
For the higher order terms, a simple, but rather lengthy calculation shows that, we have, 



d_ 

at 



r 2 fcO w , r (r,t) - 2rfcO w (r,i) - r 2 kl (l \ t (r, i)l = ^ W 



where <SW is a 'source', that depends on the lower order perturbations. We then have, 



A J = - lim 



S^+S^+S 



(4) 



dt 



(111) 



(112) 



One can show that the net contribution from vanishes. Then, the leading contribution to A J comes from . 
From Einstein's equations we find that S^ 3 ' is bilinear in the I — 2,m — ±2, order L 2 and order JL 2 contributions. 
Writing these perturbations in the Regge - Wheeler gauge, and using the expansions of the previous Section, a 
straightforward computation shows that, 



AJ=^- lim 

47T r — *-oo 



9x2 ,1 



dtp 



2,-2 



dt r2 '~ 2 dt dt 
Using now the results of the numerical evolution, we finally find, 



"X2,2 



3X2,- 



2,2 



di\) 



2,2 



Of 



-X2-2 



dt 



(113) 



AJ 
J 



= 0.0023 — 

M 



(114) 



This result implies that less than 0.3% of the total angular momentum will radiated even in the extreme case 
L ~ M. It is in good agreement with the corresponding result found in [0 for the "pseudo-inspiral" case. 
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VIII. FINAL COMMENTS AND CONCLUSIONS 



In this paper we have numerically evolved a particular set of initial data, and obtained quantitative expressions for 
the radiated energy and angular momentum, for a system of two equal mass rotating black holes, in the close limit. 
If we analyze in more detail the expression for the radiated energy, we find that a large contribution comes from the 
order J 2 amplitude. This is probably due to the fact that the initial data for rotating black holes obtained with the 
Bowen and York prescription, even for a single black hole fL5|] , contains a fair amount of "spurious" radiation that 
must be radiated away before the system can settle to the final Kerr form. One might then expect that they would be 
absent, or make a much less significant contribution if the evolution started with initial data closer to that of a Kerr 
black hole. On the other hand, it seems reasonable to assume that the contributions of the terms of order L 2 , and 
JL 2 should be much less dependent on the form of the initial data, since for order L 2 they are independent of spin, 
and for order JL 2 they require the order J data, and this is "exactly" of the Kerr form. It is precisely these terms 
that provide the leading contribution to the radiated angular momentum. This contribution is such as to decrease 
the angular momentum of the system. The radiation may therefore be considered as providing a "torque" somewhat 
in the manner envisioned in Ml. 



ACKNOWLEDGMENTS 



We are grateful to Manuel Tiglio for his help in the early stages of this work and helpful criticism and comments. 
We are also grateful to Jorge Pullin and Gaurav Khanna for their comments. This work was supported in part by 
grants of the National University of Cordoba, CONICET, Fundacion Antorchas, Agenda Cordoba Ciencia and grant 
NSF-PHY-9800973. R.J.G. is a research member of CONICET (Argentina). A.E.D. is a CONICET postdoctoral 
fellow. 



[1] See, e.g., E. E. Flanagan and S.A.Hughes, Phys.Rev. D 57 4535 (1998), and references found in this paper. 

[2] See, e.g., L. Blanchet, in Gravitational Waves, Proceedings of the Como School on Gravitational Waves in Astrophysics, I. 

Ciufolini, V. Gorini, U. Moschella and P. Fre, editors (Institute of Physics Publishing, London, 2001) 
[3] R. H. Price and J. Pullin, Phys. Rev. Lett. 72, 3297 (1994) 
[4] R. H. Price and J. T. Whelan, Phys. Rev. Letters 87 231101 (2001) 

[5] A related problem, but with spins pointing in opposite directions, has been analyzed in G. Khanna, Phys. Rev. D63 124007 
(2001) 

[6] S. Brandt and B. Bruegmann, Phys. Rev. Lett., 78, 3606 (1997) 
[7] J. M. Bowen, J. W. York, Phys. Rev. D21, 2047 (1980). 
[8] T. Regge and J. A. Wheeler, Phys. Rev. 108,1063 (1957). 
[9] F. J. Zerilli, Phys. Rev. Lett., 24, 737 (1970) 
[10] R. J. Gleiser, C. O. Nicasio, R. H. Price and J. Pullin, Phys. Rep. 325, 41 (2000). 

[11] G. Khanna, J. Baker, R. J. Gleiser, P. Laguna, C. Nicasio, H-P. Nollert, R. Price and J. Pullin, Phys. Rev. Lett. 83, 3581 
(1999). 

[12] G. Khanna, R. J. Gleiser, R. Price and J. Pullin, New Jour. Phys. 2, 3 (2000). 
[13] S. Teukolsky, Astrophys. J. 185, 635 (1973) 

[14] M. Campanelli and C. O. Lousto, Phys. Rev. D 59, 124022 (1999) 

[15] R. J. Gleiser, C. O. Nicasio, R. H. Price and J. Pullin, Phys. Rev. D 57, 3401 (1998). 

[16] D. Brill and R. Lindquist, Phys. Rev. 131, 471 (1963) 

[17] C. O. Nicasio, R. J. Gleiser, R. H. Price, J. Pullin, Phys. Rev. D 59 044024 (1999) 
[18] V. Moncrief, Ann. Phys. (NY.) 88, 323 (1974) 

[19] R. J. Gleiser, C. O. Nicasio, R. H. Price and J. Pullin, Phys. Rev. Lett. 77, 4483 (1996). 
[20] C. Cunningham, R. Price and V. Moncrief, Astrophys. J. 236 674 (1980) 
[21] K. S. Thorne, Rev. Mod. Phys. 52, 299 (1980) 



18 



